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Abstract: For points in SUSY parameter space where the sneutrino is hghter than the 
lightest chargino and next-to-hghtest neutrahno, its direct mass determination from sneu- 
trino pair production process at e^e~ colhder is impossible since it decays invisibly. In such 
a scenario the sneutrino can be discovered and its mass determined from measurements 
of two-body decays of charginos produced in pairs at the ILC. Using the event generator 
WHIZARD we study the prospects of measuring sneutrino properties in a realistic ILC 
environment. In our analysis we include beamstrahlung, initial state radiation, a complete 
account of reducible backgrounds from SM and SUSY processes, and a complete matrix- 
element calculation of the SUSY signal which encompasses all irreducible background and 
interference contributions. We also simulate photon induced background processes using 
exact matrix elements. Radiation effects and the cuts to reduce background strongly mod- 
ify the edges of the lepton energy spectra from which the sneutrino and chargino mass are 
determined. We discuss possible approaches to measure the sneutrino mass with optimal 
precision. 
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1. Introduction 

Supersymmetry (SUSY) is one of the best-motivated extensions of the Standard Model 
(SM): it stabilizes the hierarchy between the electroweak (EW) scale and the Planck scale 
and naturally explains EW symmetry breaking by a radiative mechanism. The naturalness 
of the scale of electroweak symmetry breaking and the Higgs mass places a rough upper 
bound on the superpartner masses of several TeV. This gives the Large Hadron Collider 
(LHC) the opportunity to discover the superpartners within the next years 0. Due to the 
LHC's hadronic environment, mostly strongly interacting sparticles will be produced with 
light sleptons and weakly interacting gauginos appearing in decay cascades of squarks and 
gluinos. Heavy weakly interacting states will be especially hard to detect at the LHC and 
their high-precision spectroscopy will only be possible at a future high-energy International 
Linear Colhder (ILC) g. 
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In this paper we focus on scenarios within the Minimal Supersymmetric Standard 
Model (MSSM), in which sneutrinos are lighter than the lightest chargino and next-to- 
lightest neutralino. Sneutrino decays to charged particles in the final state are then of 
higher order and therefore strongly suppressed (e.g. a 4-body decay is suppressed by an 
off-shell intermediate propagator and two additional powers of the weak coupling con- 
stant). As a result, sneutrinos decay completely invisibly into a neutrino and the lightest 
supersymmetric particle (LSP, the lightest neutralino in the considered case). A sneutrino 
appearing in a cascade at the LHC is hence lost and cannot be reconstructed. An ILC 
threshold scan is precluded for the very same reason. The only possibility to access the 
sneutrino mass in such a case is to select a well-reconstructible process where the sneu- 
trino is exchanged in a t-channel or shows up as a part of a cascade decay. The precise 
determination of kinematic distributions gives then access to the sneutrino mass. This idea 
has been proposed in Ref. Q, where it was argued that background effects are sufficiently 
under control, so a precise mass determination is possible. 

Here we study the prospects for determining sneutrino properties in a realistic ILC 
environment. We include beamstrahlung, initial-state radiation, a complete account of 
reducible backgrounds from SM and SUSY processes, and a complete matrix-element cal- 
culation of the SUSY signal which encompasses all irreducible background and interference 
contributions. Radiation, background and interference do have considerable effect on spec- 
trum shapes and edges, and we discuss the possible approaches to nevertheless determine 
the sneutrino mass with optimal precision.^ An interesting question is how well this can be 
done in a scenario with nearly degenerate charginos and sleptons, where the correspond- 
ing decay leptons are fairly soft and are in danger of being swamped by photon-induced 
background. We restrict ourselves to areas in SUSY parameter space where charginos are 
within reach of a 500 GeV ILC. 

As a working point we take the SPSla' scenario Q which has been widely in use for 
phenomenological studies. It is derived from the SPSla point Q by lowering the universal 
scalar mass parameter in order to satisfy the cold dark matter constraint. Such a shift 
leads to lower sfermion masses and, as a result, predicts sneutrinos that decay invisibly. 
The SPSla' satisfies all experimental constraints from precision data and cosmology. 

The paper is organized as follows. In Sec. § we discuss in detail the signal and back- 
ground processes, especially investigating whether they are experimentally distinguishable 
or not. Effects coming from the inclusion of initial state radiation (ISR) as well as beam- 
strahlung are discussed. The cut-based strategy to enhance the signal-to-background ratio 
is developed in Sec. ^, where we discuss the influence of the cut procedure on the quality 
of the measurement of several observables. The crucial observation here is the distortion 
of the lepton energy spectra due to ISR and photon-induced processes as decay leptons 
are quite soft in almost degenerate scenarios like SPSla'. All analyses in this Section have 



been done for a 500 GeV center-of-mass (CM) energy apart from Sec. 3.3, where important 
differences for an 800 GeV ILC are discussed. In Sec. ^, we present different approaches to 
the sneutrino mass determination in our framework. Finally, we conclude as well as give 



""^Preliminary results have been presented in j^. 
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an outlook. 

For reader's convenience, Appendix A recapitulates masses and branching ratios (BRs) 
of super symmetric particles for the SPSla' parameter point that are important for our 
analysis. Appendix B is devoted to the comparison of the effective photon approximation 
(EPA) with the exact matrix element method employed in our calculations. 



2. Signal and background processes 

After SUSY will have been found and identified at the LHC, the goal of the ILC will be the 
precision determination of the masses, quantum numbers and coupling constants in order 
to get access to the high-energy theory and possibly the SUSY breaking mechanism Q. 
Especially challenging are measurements of sneutrino parameters in scenarios in which they 
decay invisibly, as in e.g. the SPSla' parameter point. The difficulty lies not only in an 
inability of exploiting the sneutrino pair production process, but also in particularly many 
background processes, both of SUSY and SM origin, contributing to the signal signature 
as discussed below. 

2.1 Sneutrino signal and signature 

One of the standard candle processes is the pair production of the lightest chargino, 

e+e- ^ xtXi ■ (2.1) 

The chargino mass can be measured either by a threshold scan or in the continuum from 
the edges of the decay spectra In the SPSla' scenario the light charginos Xi decay 
predominantly to t^u-j- with a branching ratio BR(x^ Tj^z^t-)=53.6% (see Table A.l in 
Appendix A), followed by — > r+x?- This cascade produces a final state similar to that of 
stau pair production. Using two tau jets in the opposite hemispheres and missing energy 
as the signature for x^Xi production, the measurement of the chargino mass has been 
simulated for the SPSla scenario, and the expected accuracy is 0.55 GeV ^J. Detailed 
simulations showed that also the stau mass in SPSla can be measured with an uncertainty 



of 0.30 GeV |1C, 11]. A similar precision for the SPSla' point can be expected, since - in 
addition - chargino decays to electron (or muon) and sneutrino can be exploited. 

The main objective of our work, however, is to measure the sneutrino mass, where 
we assume the first and second generation sneutrinos {ue and u^) to be mass degenerate.^ 
This can be done by investigating the chargino decay modes Xi ~^ i^e(~^, ^ = e, /x, in 
chargino pair production process 0]. With a BR(x5'' — > U££~^)=13% a large sample of 
events can be expected in this channel. The two-body chargino decay leads to a uniform 
decay lepton energy spectrum with edges determined by the chargino and sneutrino masses, 
which provides a method for their experimental determination. 

In order to avoid large backgrounds from same-flavor lepton and slepton production, 
it is advantageous to select events with leptons of different flavor, i.e. one of the charginos 



^Note, however, that the fit presented in Sec. ^ can also be apphed when sneutrino masses of the first 
two generations are different. 
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decays to an electron, the other to a muon. So, we search for (semi-)exclusive final states: 



e+e xtXi z^eC z^^/z+ e ii^VeV^XiXi- (2-2) 

There are also other channels that contribute to the (semi-)exclusive final state as in 
Eq. Q, 

e+e- e~ ^x+ue X? X? , (2.3) 

with the intermediate state X that includes difi'erent production processes as well as in- 
terference terms, e.g. multi-peripheral diagrams, chargino decays into neutralino and W , 
smuon Jlfjlj and selectron efej pair production, and single-resonant SM di-boson pro- 
duction. In the following, we will distinguish between the "proper" signal, where the final 
state products come from chargino pairs decaying to electron and muon sneutrinos as in 
Eq. ( |2.2| ), and the signal including irreducible background, Eq. If not mentioned 

otherwise, numbers and cross sections labeled "signal" always include these irreducible 
background processes. The mixed chargino decay modes, where one chargino decays as 
xf — > i^ei~^ and the other Xi T^Ur followed by the cascade of r and leptonic r decays, 
will have additional neutrinos in the final state. Such processes are considered as reducible 
background. All reducible SUSY backgrounds that contain additional neutrinos in the final 



state will be treated in Sec. 2.2. 



At the ILC, for the SUSY spectrum of SPSla', the Born cross section for Xi P^-ir 
production without initial-state radiation or beamstrahlung is 173.56 fb at a CM energy 
of 500 GeV. The cross section reaches its maximum of ~200 fb near y/s = 600 GeV, and 
then decreases to 181.30 fb at 800 GeV, and falls further to 141 fb at 1 TeV, see Fig. | 



(tree level result shown together with the 1-loop corrections calculated in [^, 12, |T^). The 
chargino partial width for the decay into sneutrino and lepton is 10.2 MeV for each of the 
first two generations, which constitutes 13.3% of the total decay width each; branching 
fractions for other decay modes are 18.5% for tan sneutrino and tau, and the dominant 
decay with 53.6% is into stau and tau-neutrino. At 500 GeV the cross section for the 
proper signal a{e'^e~ xfxi) >^ BR{xi — ^ ^ BR{xt ~^ ^/^Z^^) is 3.06 fb, 

whereas for the six-fermion final state, Eq. (|2.3| ), it is 4.69 fb. This shows that there 
are considerable off-shell and interference effects from other non-resonant SUSY processes 
contributing to the same six-fermion final state (cf. also |jl^). These effects are additionally 
smeared when ISR and beamstrahlung are taken into account; we then obtain 2.50 fb for 
the proper signal (i.e. when forcing the final state to come from chargino pairs decaying 
to sneutrinos) and 3.94 fb for the signal. As expected, both values are slightly lower when 
ISR and beamstrahlung are taken into account, as the emitted photons drive the process 
to slightly lower values of the effective CM energy and hence lower chargino cross section. 
For a CM energy of 800 GeV, the respective numbers including initial state radiation and 
beamstrahlung are 6.60 fb for the signal and 3.23 fb for the proper signal. Here, the 
discrepancy between full matrix elements and narrow-width approximation is even more 
severe since more SUSY processes can contribute to the exclusive final state. 

Experimentally, the signature for the signal process is one electron and one anti-muon 
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Figure 1: Total cross section for chargino pair production for the point SPSla', at LO (medium 
grey, red). The NLO curves show a fixed order approach (black, blue) and with photon resummation 
(light grey, green) (from [|, |l2|). 

and missing energy,^ 

e+e" e"/i+ + ^ . (2.4) 

Such final states with two visible leptons and missing energy, however, are generated by 
many other processes, both with supersymmetric and Standard Model particles in the 
intermediate states. Many of these processes should be suppressed both by phase space 
and higher orders in the electroweak coupling constant. On the other hand, photon-induced 
processes - enhanced by large collinear logarithms - give by their sheer cross sections the 
most severe background at the ILC. In the following, we will classify the reducible and 
irreducible background processes from SM as well as SUSY productions. 

2.2 Reducible SUSY background 

Contamination from other SUSY processes comes from any process with the generic sig- 
nature 

e~^e~ any SUSY particles Xi Xi ^ /"^ + ni' , (2-5) 

with n > 4. Especially severe contaminations come from processes that lead to two lepton- 
ically decaying r's. They mainly originate from the production of mixed neutralino pairs 
Xj; stau pairs as well as chargino pairs. The latter then undergo the subsequent decays 
to stau and LSP or tau and sneutrino.^ 

In listing the contributing processes, there is a tension between describing them signa- 
ture-driven or considering them according to their exclusive production mechanisms. Since 
there are big off-shell and interference effects among several contributing weak amplitudes, 
as was shown above from the cross section considerations, it is difficult and dangerous 
to split processes completely into their corresponding production channels. (Also at the 

^Equally well, one can consider e'^ fi~ + as a final signature or use both channels. 

*Since the proposed ILC vertex detector has a granularity that is a factor 3 finer than assumed in 
Ref. one could in principle use displaced tau decay vertices and veto against explicit tau decays. This 
is a costly procedure and is to our knowledge not used in any of the ongoing experimental ILC studies. 
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Process 


500 GeV 


cj [fb], presel. 


cj™t [fb] 


Signal 


ee ^ e/xPef/.X?X? 


3.940( 8) 


1.639(3) 


SUSY r Bkgd. 


ee — > Try?Y? — > e/iY?y?4zv 


4.107( 7) 


0.978(2) 


SUSY Tu Bkgd. 


ee — > Try?y?2z/ — > e/iy?y?6i^ 


3.245(10) 


0.818(3) 


SUSY re Bkgd. 


ee — > eri/pZ/xy?y? ^ euy?y?4z^ 


3.691( 9) 


1.102(8) 


SUSY Bkgd. 


ee ^ UTfnh'^Y^Y^ — ^ ei/'y?Y?4^' 


2.617(10) 


0.966(8) 


SM WW Bked 


— >■ pulu 


152 42f25) 


736f2) 


SM er Bkgd. 


ee — > er2i^ — > e/i4z^ 


26.522(12) 


0.317(1) 


SM /XT Bkgd. 


ee — > iit2i' — > ejiAu 


15.569(54) 


0.174(1) 


SM 1/ Bkgd. 


ee — > e^Au 


0.145( 1) 


0.016(3) 


SM T Bkgd. 


ee — >■ TT — > e/i4z/ 


32.679(98) 


< 0.001 


SM ri/ Bkgd. 


ee ^ rr2i/ ^ e//6z/ 


3.852(10) 


0.335(9) 


SM 7 ^ T Bkgd. 


7*7* — > TT — > e/i2z^ 


21392(70) 


0.273(2) 


SM 7 ^ c Bkgd. 


7*7* — > cc — > e^jj2v 


1089( 4) 


< 0.001 


SM 7 ^ 1^ Bkgd. 


WW efi2u 


1.094( 6) 


0.079(1) 


SM 7 TVr Bkgd. 


7*7* tt2v — > e/x8z^ 


0.077( 1) 


< 0.001 


SM 7 ^ Bkgd. 


7*7* — > (e,/i)r2i/ ^ e/i4z^ 


0.404( 2) 


0.055(2) 



Table 1: Cross sections for all signal and background processes for an ILC energy of 500 GeV. ISR 
and beamstrahlung are always included. Note that the final states e always means electron and 
fi always anti-muon. For more details about the processes confer the text. The "presel." column 
always includes a 5° cut for the final electron to cut out coUinear regions, and for the 7-induced 
processes a 1° cut for particles vanishing in the beampipe. The last column shows cross sections 
after the cuts discussed in Sec. |[ In parentheses are the WHIZARD integration errors. 

ILC, SUSY processes appear to have a generic cascade chain structure, but due to their 
production being electroweak, a lot more interference among different chains is possible 
and actually sizable). On the other hand, an experimental description which takes into 
account only visible particles and missing energy is not much discriminating. We here 
pursue the convention to classify processes as SUSY or SM according to whether they 
come from SUSY or SM production, which is apparent from the LSP's in the final state. 

In the following we classify signal and background according to the production mech- 
anism if it is clearly identifiable (especially for the SM processes), or specifically for the 
SUSY processes we use the SM particles that appear together with the LSP in the final 
state (before a final leptonic t decay). 

In our analyses we consider as (reducible) SUSY background all processes that produce 
final states with e~/x~*", two lightest neutralinos (LSP) and four or six neutrinos. For the 
parameter point SPSla', the sum of these SUSY processes gives event numbers for the 
signal signature that are bigger than the signal typically by a factor three to four. 

We distinguish the following SUSY background processes in our study, according to 
the number of intermediate state t's and Ur^s, because t decays are technically treated 
differently (see below. Sec. \2.4( }. Furthermore, the t's are considerably long-lived such that 
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there is no interference with non-resonant diagrams, and the processes can be well split 
here. The largest contribution comes from two LSP's together with a tau pair, giving rise 
to the two tagged leptons efi, four neutrinos and two LSP's: e~^e~ — > X — > t~^t~XiXi ~^ 
^^e~ ^'^v^Vei^ri^T- The main contribution to this final state - which we call SUSY r 
background - comes from on-shell stau pair production t^^t" , as well as from neutralino 
pair X1X2 production in the intermediate state X. This process receives also smaller contri- 
butions from other intermediate X: neutralino pairs XiX^ (especially XXX3), production of 
heavy chargino X1X2 (on-shell only at 800 GeV), and even Higgstrahlung or SM di-boson 
production with (partial) decay into SUSY particles. 

Since we are exclusively looking into positively charged muons and negatively charged 
electrons, there are processes with two different intermediate states containing a single r: 
e+e~ X ^ ^i^u^T'VrXAXi ^ Xie"M"^i^^z^eZ^TZ^T, and e+e" X ^ e~ UeT^ i^rXiXi ^ 
^y^e' ^^v^Ve^T^T- For both of these processes, the main contributing intermediate states 
are chargino pairs XiXi with one of the charginos decaying directly to electron or muon and 
the other decaying to r. The other possible double-resonant contributions are, for example, 
production of neutralino pairs XiX3,45 stau pairs , and smuon Jl^JlJ or selectron pairs 
efe~, respectively. Single-resonant contributions include SM W pairs, production of heavy 
charginos X1X2 ' supersymmetrized version of vector boson fusion processes (VBF) (e.g. 
replacing the VBF final state e'^W'u by e^Xi^u) and tri-boson production. Because of 
the electron in the final state of the second process, there are many more VBF processes 
and peripheral t-channel topologies. 

The most complex process - the SUSY ti/ background - is the one which contains in 
addition to the SUSY r background also two neutrinos in the intermediate state: e~^e~ — > 
X — > T~ Ui^iXiXi ~^ Xa)^^~ ^^^^^^^el^Tl^T^^^i■, i = e,n,T, which yields 60,000 Feynman 
diagrams before the final r decays. The leading contribution is again due to production 
of chargino pairs xfxi and the following decays to r leptons. The other possibility is 
production of neutralino pair X1X2 processes similar to SUSY r background with 
additional bremsstrahlung of electromagnetic and weak gauge bosons. By far the highest 
number of diagrams is for a heavy neutralino pair production with their subsequent decays 
(it becomes important at 800 GeV, where heavy neutralinos can be produced on-shell). 

For the photon-induced SUSY processes, the largest contribution comes from stau pairs 
with a cross section of 0.028 fb at a 500 GeV ILC. Since it can be easily cut out, we do not 
consider any 7-induced SUSY background from now on. 

2.3 Standard Model background processes 

The SM backgrounds leading to the same signature of a different-flavor opposite-sign lepton 
pair and missing energy (carried away by neutrinos in the SM), mainly come from WW 
pairs, single W production and r+r" pairs. Leptonic decays of the W^s and r's lead in all 



these cases to the signal signature, Eq. (2.4). 

Pair-produced and leptonically decaying VK's are the most severe background with 
a cross section of roughly 200 fb, which is an order of magnitude larger than all SUSY 
processes. Apart from cases where the W^s decay directly to an electron and an anti- 
muon, it includes also processes in which one of the Vl^'s or both decay to tau(s), which 
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Process 


800 GeV 


cj [fb], presel. 


CT^"' [fb] 


Signal 


ee enUeUf^XiXi 


6.595(17) 


1.603(20) 


SUSY r Bkgd. 


ee — > TTy?Y? — > e/iy?Y?4z^ 


3.007( 6) 


0.731( 1) 


SUSY TV Bkgd. 


ee — > TTi^,- i/,y?y? — > e/iy?y?6z^ 


4.324(16) 


1.147( 4) 


SUSY re Bkgd. 


ee — > eTUpiy-rY^Yi 


5.458( 9) 


1.098( 6) 


SUSY T^l Bkgd. 


ee — > WTz/,,j/^y?y? 


3.639( 9) 


0.974( 8) 


SM WW Bkffd 




140 49f26) 


338f 5) 


SM er Bkgd. 


ee er2z^ — > e/i4i^ 


24.330(13) 


0.193( 1) 


SM /XT Bkgd. 


ee — > /ir2z^ — > e/x4z^ 


8.245( 9) 


0.084( 1) 


SM Bkgd. 


ee — > euVf.i'ni'ii'i 


0.214( 1) 


0.060( 4) 


SM r Bkgd. 


ee — >■ TT — > ejJLAv 


13.74 (18) 


0.005( 0) 


SM rzy^ Bkgd. 


ee tt2v e^Qv 


2.981( 9) 


0.670( 3) 


SM 7 ^ r Bkgd. 


7*7* — > rr ^ e~ ^'^2u 


28076(168) 


3.088(18) 


SM 7 ^ c Bkgd. 


7*7* — > cc ^ e~ ^'^ jj2u 


1568( 6) 


< 0.001 


SM 7 ^ Bkgd. 


J* J* ^ ^p^/ ^ e/i2z^ 


3.458(21) 


0.620( 4) 


SM 7 ^ TVr Bkgd. 


7*7* ^ rr2i/ ^ e;u8z^ 


0.116( 1) 


0.0410( 2) 


SM 7 ^ £r Bkgd. 


7*7* — > e~T'''2i/ — > e^/i^4i/ 


1.256( 6) 


0.356( 1) 



Table 2: The same as in Table but for an ILC energy of 800 GeV. 



then decay leptonically yielding our signature. Since the lepton energy spectra from the 
VF's peak at half their masses, we will later veto against very hard leptons. 

As a second severe SM background, we have leptonically decaying r pairs. Their cross 
section is 32.7 fb, being an order of magnitude larger than the SUSY signal, and still a 
factor three larger than all SUSY processes. Since the leptons originating from tau decays 
have a genuine back-to-back structure, this enables one to identify the tau background with 
a thrust-like variable. 

At the ILC there are also photon-induced SM processes that constitute a severe back- 



ground |llj, |15|, |1^. This background is usually two to four orders of magnitude larger 
than the corresponding signal one is interested in. The background stems from collinear 
photon radiation off the incoming electrons and positrons, which do not get a large kick 
and vanish in the beampipe. The collinear photons then trigger very much the same SM 
background processes, but compared to the genuine electroproduction the cross sections - 
although reduced by the square of the electromagnetic coupling constant - are enhanced 
by large collinear logarithms. 

Although the designated electromagnetic calorimeters at the ILC will have a spectac- 
ular resolution in the extreme forward and backward directions, we demand the detected 
electron in all cases to have a polar angle 9e > 5° with respect to the beam axis. This cuts 
out the collinear singular region, where perturbation theory is no longer reliable. For the 
photon induced processes, particles with a polar angle < 1° with respect to the beam di- 
rection will be treated as lost in the beampipe. Therefore, we preselect events by requiring 
6e > 5°, while for other charged particles considered to be lost in the beampipe - jets and 
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additional leptons - we demand ^ < 1°. 

While the total cross section for photon induced processes is huge (about 12 nb in 
total, about 900 pb for charm, 300 pb for tau, and 90 fb for W pair production), it greatly 
reduces with the decay branching fractions to the required final state folded in, and after 
preselection. Nevertheless, these processes are still very large, amounting to roughly 23 pb 
at a 500 GeV ILC. By far the most dominating background is from photon- induced tau 
pairs decaying to the tagged leptons with a cross section of 21.4 pb. Their broad spectrum 
at low energies together with soft leptonic decays distorts the shapes of the low-energy edges 
of the electron and muon distributions that are important for the chargino and sneutrino 
mass determination. 

Photon- induced charm pairs have even larger total cross section than the r's, but 
vetoing against jet activity in the central detector from their semi-leptonic decays cuts this 
down to a value of 1.1 pb, which is still a factor 250 larger than the signal. Leptonically 
decaying Ws contribute another 1.1 fb. 

The cross sections for the signal, the SUSY and SM backgrounds after preselection 
for ^/s = 500 GeV ILC are collected in Table |^. The right-most column in Table |^ shows 
the change in cross sections after applying suitably chosen cuts which will be discussed in 
detail in Sec. y. The corresponding values for y/s = 800 GeV are given in Table || 

2.4 Event simulation 

For the simulation of the SUSY signal and SUSY/SM background processes, we take the 
multi-purpose event generator WHIZARD [o], which is especially suited for beyond the SM 



applications and well-established for SUSY simulations [18|. It allows for the usage of 
full matrix elements for exclusive final state particles, and automatically generates all 
contributing intermediate (on- and off-shell) states. Hence, intermediate resonances as 
well as interferences from off-shell continua are considered on equal footing. For ILC, the 
inclusion of off-shell states in full matrix elements is mandatory, especially when cuts have 
to be taken into account |p^]. 

Since the process considered here (as most electroweak SUSY processes and SUSY 
decay cascades) are dominated by tau leptons, it is crucial to simulate leptonic decays of 
taus. Since taus are very narrow resonances, a multi-channel adaptive phase space routine 
has difficulties finding the corresponding poles. In that case, a narrow-width approximation 



for the tau decays is appropriate [19|. We simulated the energy and angular distributions of 
the leptonically decaying r's and the semi-leptonically decaying charm quarks by extending 
the WHIZARD program accordingly. 

As mentioned above, the most severe background at a high-energy lepton collider comes 
from photon- induced processes. Such contributions are commonly accounted for by using 
the equivalent photon approximation (EPA) |2^, which approximates the collinear ra- 
diation by a structure function approach of on-shell photons inside the electrons. However, 
the EPA tweaks the kinematics and results in deviations in total cross sections as well as 
in the shapes of differential distributions. This is particularly important when kinematic 
cuts are used to select the desired signal. Taking the photon-induced tau-pair produc- 
tion process as an example, in Appendix B we compare EPA with the exact treatment 
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exposing the difference in the lepton energy spectra before and after cuts. In WHIZARD, the 
approximation of on-shell intermediate photons is avoided by using exact matrix elements 
for cohinearly radiating electrons that vanish in the beam pipe. However, it is necessary 
to run WHIZARD in a higher (quadruple) floating point precision for taking properly into 
account the large p± logarithms of the radiating electrons. 

Since all distributions at the ILC are heavily dependent on photon initial state radiation 
as well as beamstrahlung effects, it is absolutely mandatory to include these effects. This 
is done using the setup inside the WHIZARD generator, which employs an ISR structure 
function that resums all leading soft- and soft-collinear logarithms and takes into account 
hard-collinear terms up to the third order [22|. The beamstrahlung is simulated online 
using the CIRCE package [^] that comes with the WHIZARD distribution. In ||l^ it has 
been shown that ISR and beamstrahlung can heavily distort the box-shaped spectra of 
energies from decay products of heavy states whose mass one wants to determine. Since 
at 500 GeV the collider energy is at the rising shoulder of the production cross section for 
charginos in the SPSla' scenario, the radiative return slightly reduces the number of signal 
events, bringing the signal process closer back to the threshold. On the other hand, the 
SM background, dominated by the WW production, is slightly enhanced by the radiation 
effects. Therefore in the following, beamstrahlung and ISR are always taken into account. 

Although we did not study effects of final state radiation (FSR), we briefly want 
to comment on its possible effects. Prom fast quasi-massless fermions like electrons and 
muons, final state collinear photon radiation could be enhanced by log {Ei/5i), where 5i is 
the energy resolution for collinear photons. Since for the near-degenerate SPSla' spectrum 
the decay leptons are rather low-energetic, we do not expect this to be a sizable effect. 



3. Cut-based strategy for enhancing the signal to background ratio 

Since the final state e~ fi^ + ^ contains only two visible particles, it might seem to provide 
only a limited amount of experimental information. Nevertheless, several experimental 
observables can be exploited to enhance the signal to background (S/B) ratio by applying 
suitable cuts on the lepton polar angle 6{i) (measured from the e~ beam direction), on en- 
ergy E{i) and transverse momentum p±{£) distributions, on azimuthal angular separation 
of outgoing leptons A(/>(e,/i), as well as on missing energy and momentum. 

3.1 Cutting out the SM and SUSY background 

The cuts to enhance the S/B ratio have been optimized for a CM energy of 500 GeV, with 
slight modification for 800 GeV (cf. Sec. 

As was discussed in Sec. |2^, after preselection the most severe background comes from 
photon-induced SM processes which are larger than the signal by a factor 10^. Enhanced 
by large collinear logarithms, these processes contribute to the final state leptons from the 
leptonic r and semi-leptonic charm decays. Hence, e and fi appear soft and along a thrust- 
like axis parallel to the beam with a large azimuthal separation. Therefore, a cut on the 
total transverse momentum \p±{e)+p±{fi)\ > 4 GeV (and also on the separate p_L(e, /i) > 2 
GeV) together with a cut on the azimuthal separation of the leptons |A(/)(e,/i)| < 150° 
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^/s, Cut 


500, lower 500, upper 


800, lower 800, upper 


e(e-) 

A</>(e-,/i+) 

P±(e~), 

|pj_(e") 


1 LieV 40 LreV 
15° 155° 
25° 165° 

-150° 150° 

2 GeV 
4 GeV 


1 GeV 60 GeV 
25° 165° 
10° 165° 

-150° 150° 

2 GeV 
4 GeV 



Table 3: Cuts used in the analysis for an ILC with 500 GeV and 800 GeV CM energy, respectively. 

provide a very powerful means to reduce the gamma-induced background by roughly two 
to three orders of magnitude. [Note, that a different strategy was followed by the LEP 
experiments, namely to determine the gamma-induced background from a control sample 
by tagging the emission of an additional hard photon [^].] In addition, the azimuthal 
correlation cut on the e, \x system eliminates the SM r background, since the r's are heavily 
boosted. 

Background from VK-pair production contributes to the higher-energy parts of lepton 
spectra, since the decay leptons shape a Jacobian peak at approximately half of the W 
mass. The considerable boost of the shifts this peak further up. This argument applies 
to both direct leptons as well as to e, [i coming from PF's going into leptonically decaying 
r's. Therefore, constraining final state leptons to the energy window 1 < i?(e, /z) < 40 
GeV reduces the W backgrounds (which is the largest SM background after cuts), by two 
orders of magnitude. 

Furthermore, we require two well-visible leptons in the central part of the detector, 
i.e. 15° < d{er) < 155° and 25° < 9{ii^) < 165°. The complete set of cuts reduces the 
SM backgrounds by four orders of magnitude, yielding a S/B ratio of order 1. All cuts are 
summarized in Tab. 0, and their effects on signal and background processes are collected 
in Tab. In the latter table cross sections are shown in units of fb. In the second column 
cross sections after preselection are shown, in the following columns when relaxing the 
corresponding cut(s). The most efficient cuts for the individual processes are indicated by 
shading the corresponding entry in grey. The last column shows the result after applying 
all cuts. Note that a cut on missing energy or missing transverse momentum does not help 
in enhancing further the SUSY to SM ratio since there is severe background from Z ^ vv 
and vector boson fusion. 

Next, we investigate how the cuts affect the SUSY processes. There are two issues here: 
one is the separation of the signal process from reducible SUSY backgrounds, the other 
is a desired enhancement of the rate of proper signal over irreducible SUSY backgrounds. 
The latter is important, since we want to apply kinematic analysis to determine sneutrino 
and chargino parameters, which mostly rely on the on-shell kinematics of the proper signal 
process (cf. Sec. ^). 

SUSY cascade decays for parameter points with an almost degenerate spectrum of 
weakly interacting states, like SPSla', generate relatively soft leptons with an energy of 
a few up to few tens of GeV. Since the leptons here originate from decays of quite heavy 
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Process 


presel. 


no 


no E 


no 9 


no 


all cuts 


Signal 


3.940 


1.715 


2.425 


1.842 


2.419 


1.639 


SUSY r 


4.107 


1.570 


1.327 


1.021 


1.536 


0.978 


SUSY Tu 


3.245 


1.545 


0.964 


0.851 


1.017 


0.818 


SUSY re 


3.691 


1.568 


1.429 


1.152 


1.557 


1.102 


SUSY r/i 


2.617 


1.332 


1.052 


0.999 


1.375 


0.966 


SUSY bkgd. 


13.66 


6.015 


4.772 


4.023 


5.485 


3.864 


SM WW 


152.42 


0.762 


40.98 


1.023 


0.926 


0.736 


SM er 


26.522 


0.372 


6.61 


0.500 


0.387 


0.317 


SM 


15.569 


0.195 


3.064 


0.202 


0.228 


0.174 


SM T 


32.679 


0.009 


< 0.001 


< 0.001 


1.353 


< 0.001 


SM TV 


3.852 


0.431 


0.885 


0.433 


0.535 


0.335 


SM 7 ^ T 


21392 


615.9 


0.273 


0.371 


23.54 


0.273 


SM 7 ^ c 


1089 


0.216 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


SM 7 ^ Vt^ 


1.094 


0.079 


0.628 


0.082 


0.091 


0.079 


SM 7 ^ £t 


0.404 


0.062 


0.114 


0.062 


0.067 


0.055 


total SM 


22714 


618.0 


52.55 


2.673 


27.13 


1.696 



Table 4: Effects of cuts on the signal and background processes at CM energy 500 GeV. All 
values are given in fb. Description of columns is given in the text; the last column corresponds to 
the cuts of Tab. S. The most important cuts for the individual processes are shaded in grey. 



states, they are radially distributed in the detector, and electron and muon are almost 
angular uncorrelated. Due to the complexity of electroweak SUSY amplitudes and their 
corresponding phase space configurations, there is in general no single well-motivated cut 
for reducible SUSY background. Interfering cascades, especially semi-resonant processes, 
with only small differences in kinematics are difficult to cut out. Nevertheless, with the set 
of our cuts these backgrounds are reduced by a factor of three to four. 

The irreducible SUSY background is also reduced, although to a lesser extent: a 40 % 
contamination from non-resonant or partly resonant contributions is brought down to a 
mere 13%. This is more pronounced for the muon distribution (Fig. |2|) than for electrons 
in the energy range left over after the cuts. Most of the irreducible background, which 
comes from multi-peripheral i-channel diagrams, is cut out by the upper energy cut of the 
e~ . The rate of the proper signal itself suffers mostly from the A(j), but only moderately. 
Nevertheless, after all cuts we have a clearly visible SUSY signal with S/B ~ 3/1. 

3.2 Cut effects on lepton energy distributions 

For determination of sneutrino and chargino masses from kinematic variables at the ILC, 
the energy distributions of the decay leptons are crucial. It is important to realize that these 
are distorted by ISR, beamstrahlung and selection cuts, as well as by contributions from 
background events that survive cuts. The aim of this section is to discuss the consequences 
of these distortions. When plotting only the proper signal (i.e. restricting to intermediate 
charginos decaying to sneutrinos), the famous on-shell box shape of the lepton spectrum 
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Figure 2: Electron (upper row) and muon (lower row) energy distributions for the proper 
signal (dark grey, blue) and full signal with irreducible SUSY background (dark+medium grey, 
blue+green). Left panels show histograms after preselection, and right panels after all cuts. In the 
muon distribution, the ratio of proper over full signal is considerably enhanced by the cuts. 
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Figure 3: Lepton energy distributions after applying the cuts of Tab. ^: SM background (light 
grey), SUSY background (medium grey, red) signal (dark grey, blue) for the electrons (left panel) 
and muons (right panel). 



is clearly visible, Fig. ||.^ The cuts applied to enhance the S/B ratio do not do much harm 
to this feature: although smeared to a certain degree (more for electrons than for muons) 
the box character is still visible. 



^In all histograms the number of events per 1 GeV bin assuming an integrated luminosity of 1 ab 



shown. 
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Figure 4: Lepton energy distributions for the signal (dark grey, blue) and reducible SUSY re 
background (medium grey, red). Upper row: electrons after preselection (left panel) and with all 
cuts (right panel). Remaining plots are for muons: after preselection (middle, left), after energy 
and angle cuts (middle, right), after energy, angle and p± cuts (lower, left), and after all cuts (lower, 
right). The low-energy edge of the box- type energy distribution of muons is no longer visible after 
the p± cut. Note different vertical scales in these plots. 



The lepton spectrum from the SM background after cuts is rather flat, the light grey 
histograms in Fig. ^, which allows one to subtract it in the analysis. The SUSY signal 
sits nicely on top of it, making a distinction between SUSY and SM quite easy. On the 
other hand, it is not possible to discriminate between SUSY background (reducible and 
irreducible) and the proper signal with any kind of cuts. Especially the similarities in the 
kinematic set-up of SUSY r and SUSY t£ + tv backgrounds make it difficult. 

To expose our point better, let us consider a specific SUSY re background process 

e+e" — > xtXi — ' e'Vlr+Vr — > e"PeM+ff»z^TZ^TXiXi , (3.1) 

where one of the charginos decays directly into e and t'e, while the final muon is the decay 
product of the tau coming from the other chargino. The same discussion applies to the 
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case with e, // flavors exchanged, and for the other r-induced SM or SUSY background 
processes. 

Figure |^ shows the lepton spectra for the full signal, Eq. ( |2.3D (dark grey, blue), and 
for the SUSY re background, Eq. ( |3.1| ), (medium grey, red); in these plots SM background 
has been removed. For the electron spectrum the SUSY re background does not impair 
the edges. On the contrary, it enhances the number of signal-like events, see the first line of 
Fig. This should not be astonishing since the SUSY re (and Tfi) background processes 
share in one hemisphere the same chargino cascade decay chain with the signal process. 
As a result, for electrons the kinematics in both cases is almost identical. 

For the muons in this particular process the situation is quite different. Here the 
spectrum at the lower edge of the box eats into the raising infrared tail of the SUSY re 
background. The middle and lower lines of Fig. |^ show the successive application of the 
cuts listed in the previous section to eliminate the SM backgrounds. After the p± and 
A(j) cuts, the lower edge of the muon spectrum is extinct; releasing the p± cut, however, 
reintroduces a large SM background (cf. Sec. 3.1). 



The complete lepton spectra with signal, remaining SUSY background as well as SM 
backgrounds are shown in Fig. |3|. The lower edge, which in the signal is due to the sneutrino 
decay, is essentially smeared out. The upper edge is more pronounced. We will discuss 



how to determine the sneutrino mass with this complication in Section |4J . 
3.3 Increasing the collider energy 

It is interesting to investigate possible benefits of increasing the collider CM energy to 800 
GeV. In this case, the SPSla' scenario exposes different phenomenological features than 
at lower CM energy. On the one hand, the chargino pair production is bigger by a factor 
of about 1.5 with respect to 500 GeV, see Fig. |^, which increases the signal statistics. At 
the same time, however, the irreducible background is also bigger, since now many new 
SUSY intermediate states are either above or close to their thresholds. The same argument 
applies to most of the reducible SUSY background, while the f pair production is slightly 
smaller. The SM background is also reduced: the WW cross section falls by 10%, Tfi by 
a factor 2, while the r pair production cross section is smaller by a factor of three. On 
the other hand, we are now on the falling shoulder of the chargino pair production cross 
section, so we experience larger effects of beamstrahlung and ISR similar to a radiative 
return effect. The radiative return also enhances the photon-induced rr by 30% and WW 
roughly by a factor 3. 

The result for all SUSY signal and background processes for a center of mass energy 
of 800 GeV before and after cuts can be found in Table ^. Note that the cuts here are 
adapted to the higher energy: the upper edge of the lepton energies was raised to 60 GeV, 
since the upper edge at this energy is at £' ~ 45 GeV. Also less restrictive polar angle cuts 
for muons have been applied, see Table ^ for details. Comparison of last columns of Tables 
|l] and H shows that increasing the ILC energy does not change cross sections after cuts by 
a large factor. They are of similar order. 

More dramatic change, however, is observed in the lepton energy distributions. With 
increased collider energy the infrared and collinear logarithms are much more pronounced. 
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Figure 5: Lepton energy distributions for electrons (left panel) and muons (right panel) after 
applying the cuts of Tab. || for = 800 GeV; SM background (light grey), SUSY background 
(medium grey, red) and the signal (dark grey, blue). 
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Figure 6: Muon energy distribution for the proper signal (dark grey, blue) and irreducible 
background (medium grey, green) at = 800 GeV; left: without cuts, right: with cuts. It is 
clearly visible that the ratio of the irreducible background to proper signal is greatly reduced. 



As a result, the low-energy tails are basically parallel for SM and SUSY, which means 
especially that the SM background from photon-induced r's is no longer flat. As can be 
seen in Fig. ^, after the cuts from Tab. ^ all spectra look quite similar. Although after cuts 
the background is reduced by a factor 10^ and the S/B ratio is still quite good ~ 1, the 
shapes of the distributions are much less favorable. The exponentially rising low energy 
structure from leptonic r decays overlap with the crucial lower edge region. Moreover, not 
only the lower, but also the upper edge is no longer clearly detectable. 

To understand the origin of stronger smearing effect for the upper edge of energy distri- 
bution, let us assume that the SM background can be measured well enough (presumably 
from continuum measurements below/far away from SUSY thresholds) and subtracted. 
The muon energy distributions for SUSY processes only are shown in Fig. ^. Before cuts 
the box-shape character is still visible although the upper edge is less pronounced than at 
500 GeV. Roughly half of the signal consists of irreducible background, which goes down to 
about 25% after applying cuts to suppress the SM background. However, at the same time 
the upper edge is killed. This is mainly due to the cut. With increased collider energy 
charginos are more boosted producing more events with back-to-back leptons. Such events 



-16- 



are removed by the A(/> cut, which however cannot be relaxed because of the necessary 
suppression of the huge background from photon induced r pair production. A similar 
effect, however not that strong, is visible in Fig. ^ at CM energy 500 GeV. 

The above discussion shows the importance of adjusting the collider energy to the 
spectrum. Increasing the ILC energy is more desirable to detailed studies of heavier SUSY 
particles. 

4. Sneutrino mass determination 

4.1 Sneutrino mass determination from total cross sections 

In the chargino production and decay process e"*"e~ xfXi ~^ e^/i+ + the dependence 
on the sneutrino mass enters both in the production and decay matrix elements. In the 
production process the electron sneutrino is exchanged in i-channel^. For SPSla' point this 
leads to moderate dependence of the cross section on this parameter, which amounts to a 
1% change of the cross section for a 1 GeV change in the sneutrino mass at ^/s = 500 GeV. 
At the same time the sneutrino mass enters the decay matrix element, since the two-body 
decay channel xt ~^ ^^£^~^ is open. For the SPSla' scenario approximately 50% of charginos 
decay to one of the sneutrinos (cf. Appendix 

Specifically, in the current analysis we concentrate on chargino decays to electron and 
muon sneutrinos. We find that for this particular final state, Eq. (|2.2D , the sensitivity 
of the production cross section (together with irreducible background) on the sneutrino 
mass amounts to ~ 25%/GeV, which is much stronger comparing to the sole production 
process of charginos. This is the consequence of strong kinematic effects in the decay am- 
plitude when the parent and daugther particles, chargino and sneutrino, are close in mass. 
Therefore looking at this particular channel gives a very good opportunity to determine 
the sneutrino mass. 

Since experimentally we cannot distinguish chargino production and decay process 
from the reducible and irreducible SUSY background, the latter also have to be included 
in the analysis. This will dilute the strong effect coming from the signal. However, since 
some of these background processes (e.g. SUSY re, Tfi and rz^) also involve chargino decays 
mediated by on-shell sneutrinos, we can still expect quite strong dependence. Another 
major source of uncertainty that has to be taken into account in determining the value of 
the cross section, apart from the sneutrino mass, comes from the error on chargino mass. 
Assuming that the chargino mass can be measured elsewhere with a precision of 1 GeV, it 
gives an uncertainty of signal and SUSY background cross section after cuts of order 0.7 fb 
(i.e. ~ 17%/GeV). This high sensitivity can again be attributed to the kinematic effects 
appearing in the chargino decay matrix element. 

A fit of the collider data to different samples of Monte Carlo data for varying (electron 
and muon) sneutrino masses can thus be used to get an estimate of the sneutrino mass. 
Figure |^ shows the number of events after all cuts for the signal (left) and for signal and 
SUSY background (right). We assume that the SM background, which is fiat, is understood 

®This dependence can also be exploited using forward-backward asymmetries, see for details. 



-17- 



^-events 
4000 



#'evonts 





Signal after cuts 




^ ^ 500 GcV 


I. 


/ £ - 1 ab"^ 


'I, 












... . 


nominal value 


^ 



174 

n,> [GeV] 



all SUSY after cuts 

= 500 GcV 
/ £ = 1 ab"' 




Figure 7: Total number of events after cuts as a function of sneutrino mass for signal only (left) 
and with SUSY backgrounds included (right). The horizontal lines give 1-sigma statistical error for 
the nominal SPSla' sneutrino mass (172.52 GeV) and for the integrated luminosity C — 1 ab^^. 



well enough and can be subtracted. Our "collider data" are the MC sample for the nominal 
SPSla' mass. The characteristic dependence on the sneutrino mass is apparent, although 
the curve becomes flatter when including the SUSY background. The SUSY background 
has been generated with a varying sneutrino mass as well. Taking into account statistical 
error and the uncertainty due to the chargino mass error, this simple method allows to 
determine the sneutrino mass range of roughly 

~ 172.50 ± 0.75 GeV . (4.1) 

No other systematic or parametric errors have been included here. In principle, generation 
of the MC samples with varying sneutrino masses demands the knowledge of considerable 
parts of the SUSY parameters. This inverse problem can be solved by fitting the experi- 
mental data in global fits to the corresponding models Nevertheless, our simplified 
method already allows us to narrow down the possible region of the sneutrino mass. A 
more spectrum-/model-independent way of sneutrino mass determination using kinematic 
methods is discussed in the next subsections. 

4.2 Sneutrino mass determination from on-shell kinematics 

Pure kinematic relations can be exploited to determine the masses of sneutrino and, even- 
tually, chargino in a model-independent way. Masses of these two superpartners can be 
determined from kinematics of on-shell chargino production followed by the on-shell 2-body 

. For the given incoming e'^e~ energy the lepton decay spectra are 
uniform between the two energies Emm and -Emax; given by: 

where = (m~ — m~) / (2m^) is the lepton energy in the chargino rest frame and f3 = 
(1 — 4m~/s)^/^ is the chargino velocity in the CM system. Inverting, we get 

^xt = ^' E +E ' "^^^="^xfV ^^^^ 



Xi 



' decay 
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The chargino polarization effects |Q, which are all included in our calcula- 

tions, affect the uniform distributions and their response to kinematical cuts. Beam- and 
bremsstrahlung also modify the shape of the energy distributions. Other sources of smear- 
ing come from off-shell and interference effects, which however are not overwhelmingly 
important since for the SPSla' point charginos and sneutrinos are very narrow. Finally, 
irreducible SUSY and SM backgrounds, as discussed in the previous section, make the 
determination of the edges of lepton energy spectra more difficult. 

A naive estimate of the position of the edges (or a corresponding fitting to the box- 
shape) could only be reliably performed if the feature is significantly distinct. Including 
all backgrounds, this is not the case. While the position of -Emax is rather well visible, 
reading off -Bmin is questionable since the low edge is essentially smeared out by the SUSY 
background. The main culprit is the irreducible SUSY background processes, where one of 
the observed leptons (e or n) comes directly from chargino decays into iu£ and the other 
from the r that is the decay product of the chargino or even the stau, i.e. xt ~^ t'^i't or 
Xi t'^Vt ^ 'T^t^rXi discussion in Sec. 3^). The energy spectrum of the tau 

decay leptons is peaked at low energy. The transverse momentum cuts, which are necessary 
to suppress the huge photon-induced SM background, tweak these low-energy distributions 
in such a way that the -Emin edge becomes completely swamped. 

If the chargino mass, however, was determined elsewhere (from e.g. the threshold 
scan in a clearly visible and distinguishable channel), the £^max can be used to derive the 
sneutrino mass with the help of 

,4.4, 

Taking the chargino mass as an input with a conservative estimate on its error, m~± = 
184 lb 1 GeV, as well as a rough estimate for the upper lepton edge, Smax = 24 it 2 GeV, 
we obtain 

mc = 173.1 ± 1.3 GeV, (4.5) 

which agrees within large error with the input value of the sneutrino mass = 172.52 
GeV. The large error is the direct consequence of the large error for -Emax- 

For the purpose of the above analysis the edge -Bmax and its error has been determined 



"by eye" |3C]. It is common when extracting the observed edges from plots, such as those 
above, to fit a function to the endpoint in order to determine both the precision and the 
accurate position. However, when taking into account ISR, beamstrahlung, background 
processes, off-shell and interference effects and their response to kinematic cuts, the form of 



fitting function is impossible to determine analytically. As stressed in Ref. [SI], even in fa- 
vorable cases applying analytic functions too readily may lead to inaccurate measurements 
and underestimated errors, since endpoints can often exhibit tails or get smeared. There- 
fore, in the next subsection we develop a different experimental method of determining 
-E-max that allows for a more precise measurement of m^. 

4.3 Sneutrino mass determination from x^ bin by bin fitting 

In order to use maximum information from energy distributions of the final state leptons. 
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we apply another method of determining the sneutrino mass at y/s = 500 GeV, where we 
perform a binned fit- This is accomphshed by using the MC true sample for SPSla' as 
the "experimental data" and generating MC control samples for varying sneutrino masses 
within the range 171 < mi> < 174 GeV in a small window around the true value. For this, 
we assume knowledge of the approximate sneutrino mass range from one of the methods 
presented in the previous two sections. Despite the fact, that this method can also be 
used in the case of non-degeneracy in the first two generation sneutrinos, we keep that 
assumption for the following analysis. In a non-degenerate case, the fitting had to be 
applied for each lepton energy distribution separately. Here, however, we combine the 
information from both histogram fits. The tau sneutrino mass is left at its nominal SPSla' 
value. For the histograms, we use a conservative value for the binning of 1 GeV. The signal 
and all SUSY backgrounds are considered for the analysis. 
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Figure 8: x^/dof distributions for the sneutrino mass determination. Left: fitting the whole 
energy spectra of both e and /i, right: fitting the bins around the upper edges. Dashed line are 
with, soUd hues without rescaling (see text for explanations) . Horizontal lines denote a significance 
level of 10%. 

Mathematically, is the best estimator whenever the uncertainties of each bin are 
Gaussian. Since this is the case (Poisson statistics from Monte Carlo with large enough 
event numbers), we use a fit the mass determination. We are aware of the fact that 
matrix element methods give the best statistical means for a particle mass determination 
(as e.g. for the top mass determination at Tevatron and LHC). However, such an analysis 
is beyond the scope of the present paper. The x'^ statistics used in our fit has the following 
form 

(iVf -iVf°°(m^))2 



iVf + iVr°(m^) 



(4.6) 



where the sum runs over bins in energy distributions of electron and muon, A^^*^^(?7t.£>) is 
the number of expected events in the i-th bin from our control sample as a function of 
the sneutrino mass, and Nj^'^ is the number of observed events in the i-th bin for the true 
SPSla' sneutrino mass; statistical errors are added in quadrature. Using Eq. (|4.6| ), we 
calculate x^ statistics for each of the control sample histograms. 

We perform two types of fits. In the first, we use all bins in the full energy range 
2 — 40 GeV, yielding 76 degrees of freedom (dof) when combining electron and muon 
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data. In the second, we use only bins in the vicinity of the upper edge, which is assumed 
to be experimentally determined - at least approximately - by e.g. a total cross section 
measurement (see Sec. 4.1). As a window, we use ±4 GeV around the expected upper 
energy edge of the decay leptons. The second method amounts to 18 dof for the 
statistics. We perform each of the two fits in two variants: using the absolute numbers of 
events per bin in the control sample, and in a rescaled version. In the latter we normalize 
the control sample to the number of true events, i.e. we define for each bin 



con, rescaled 



n: 



con tot 



ATcon 



(4.7) 



where N^^^ and N^°^ are the total number of events in the true and control samples, 
respectively. Since the normalization to the true sample imposes an additional constraint, 
it decreases effectively the number of dof by one. In this case the is sensitive only to the 
shape of energy distributions which does not rely much on the knowledge of the rest of the 
SUSY spectrum. Therefore, the error for the sneutrino mass determination is dominated 
by statistics, and does not suffer too much from parametric uncertainties. 



xVdof 




xVdof 




171 171.5 172 172.5 173 173.5 174 

m,> [GeV] 

Figure 9: x^/dof distributions for the sneutrino mass determination. Fits for the complete energy 
spectra (solid) and in the window around the e and /x edges (dashed lines). Left panel for rescaled 
number of events in the bins, right panel without rescaUng (see text for explanations). Horizontal 
lines denote significance level of 10% for 76 (solid, full range) and 18 (dashed, restricted window) 
dof, respectively. 



The plots in Fig. ^ for the x^/dof calculated according to Eq. ( ^I^ ) compare the two 
methods with and without rescaling the number of events in each bin. On the left hand 
side, the fits over the whole energy range are shown, whereas the right plot shows the result 
of restricting the fit to the window around the edge. The horizontal lines denote a 10% 
significance level for the x^/dof. In Fig. |9| we illustrate the difference between x^/dof in 
the full energy range versus the restricted energy window fits with (left panel) and without 
(right panel) rescaling. 

As can be seen from both Figures § and|^, the distributions show significant statisti- 
cal fiuctuations, especially around their minima. These arise due to statistical fluctuations 
in the generation of the MC control samples, as well as those in the true MC sample, where 
the latter would correspond to the true experimental statistical uncertainty. 
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In order to assess the possible accuracy in the sneutrino mass determination, we test 
the hypothesis whether the control sample for a given mj; coincides with the MC true ( "ex- 
perimental") data. We define the compatibility region for the sneutrino mass as the range 
where the null hypothesis cannot be rejected at the significance level of 10%. Note, that this 
approach does not directly translate to a confidence level (CL) interval for the sneutrino 
mass, but rather gives a conservative estimate. Since we did not include any systematic 
and parametric uncertainties, we consider such a conservative estimate a reasonable alter- 
native to the usual procedure of taking the respective region around the minimum of the 
distribution. 



method — > 


full range 
rescaled 


full range 
not rescaled 


edge 
rescaled 


edge 
not rescaled 


limit I 


lower 


171.9 


172.1 


172.1 


172.2 


upper 


173.1 


173.0 


172.9 


172.8 



Table 5: Upper and lower limits of the compatibility range for the sneutrino mass (in GeV) 
corresponding to the significance level of 10% of testing the null hypothesis using different methods 
for fits. The true value for the sneutrino mass is m^r — 172.52 GeV. 

The results of fits using different approaches discussed above are collected in Table |5|. 
In general, the resulting accuracy is of the order of 0.5 GeV. The best precision is obtained 
for edge fitting without rescaling samples. Edge fitting quite generically gives better results. 
Although, when using the full range of all 76 bins in the fit, we utilize more information 
about the distribution, a lot of it is irrelevant for sneutrino mass determination, hence 
giving a somewhat broader minimum of the function. As mentioned before, since we 
loose some information concerning the total cross section (or, on the contrary, dependence 
on the SUSY model) when rescaling data, therefore - as expected - in this case the fits 
give lower precision. On the other hand, rescaling minimizes effects which can result from 
parametric uncertainties on the SUSY model. 

In summary, we studied three different methods of sneutrino mass determination, in- 
cluding all SM and SUSY backgrounds. A first estimate, assuming the rest of the spec- 
trum is known, can be obtained from the total number of SUSY events (signal and all 
backgrounds) after the application of cuts. 

Alternatively, we can use purely kinematic information from the two subsequent two- 
body decays of the chargino and the sneutrino; in that case, the error is twice as large 
as the one obtained from the total cross section, but does not depend as strongly on the 
specific SUSY model. The main difficulty of this approach is defining the position of the 
upper edge, which is smeared by cut effects, ISR and beamstrahlung. 

The most sophisticated method that we developed employs a fit of the "experimental" 
lepton energy spectra by template distributions generated with varying sneutrino mass. 
The method allows one to include or exclude additional knowledge on the SUSY model. 
Here, the precision of the sneutrino mass determination can be reduced to ±0.4 (±0.3) GeV, 
when performing edge fits without (with) additional information on the SUSY spectrum. 
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respectively. We used the simplest version of this method; a full analysis would require 
matrix element reweighting for a correct treatment of the statistical properties of "true 
data" and "control samples" with varied sneutrino masses. On top of that, a treatment 
of systematic uncertainties (like detector effects) is mandatory to get a decisive answer 
on the final precision of the sneutrino mass determination. This should be a part of an 
accompanying experimental study. 



5. Conclusions and Outlook 



In this paper we have investigated the problem of the discovery of the invisibly decaying 
sneutrino and its mass determination in a realistic ILC environment. We considered the 
SPSla' scenario in which sneutrinos are lighter than the lightest chargino and next-to- 
lightest neutralino. Since decay modes with charged particles in the final state are of higher 
order, and thus strongly suppressed, sneutrinos decay invisibly via a tree-level process into 
the lightest neutralino and neutrino. As a result, sneutrino masses cannot be measured 
by a threshold scan. Therefore we have analyzed the opportunity to measure sneutrino 
masses from the kinematics of the chargino pair production at the ILC followed by two-body 
chargino decays xf — > £^i)^*^ (with i = e, /x). 

The two-body chargino decays generate sharp edges in the lepton energy spectra which 
depend on chargino and sneutrino masses. The edges, however, get smeared by ISR and 
beamstrahlung. Moreover, the signal process is overwhelmed by a plethora of background 
processes generated by both SUSY and SM mechanisms. Since the selection procedure 
required to enhance the signal/background ratio distorts the lepton spectra from the signal 
and from various background processes in different manners, particular attention has been 
paid to the pollution coming from these effects. 

We have considered all SUSY and SM background processes that could be of relevance 
to the analysis. The main observation is that the leptons coming from the leptonic decays 
of the final r's (copiously produced by SUSY and SM processes) practically wash out the 
lepton low-energy edge preventing precise determination of both the chargino and sneutrino 
masses. If the chargino mass, however, is predetermined from a different observable, the 
high-energy edge alone is sufficient to measure the sneutrino mass with an accuracy of 
order a GeV. Using a fit of control samples generated for varying sneutrino mass to the 
"experimental" one at the nominal value we find that the error can be further reduced by 
a factor 4. 

The calculation has been performed at tree-level, demonstrating the capability of 
WHIZARD to generate exclusive final states using full matrix elements for production and 
decay, including off-shell and interference effects. The expected high experimental preci- 
sion at the ILC, however, calls for more precise theoretical methods. Recently the chargino 
pair production process at one-loop has been incorporated into WHIZARD |12|. One-loop 
calculations of the chargino decay process are also available |32], so a full NLO calculation 
of the signal process, which has to be incorporated in the simulation, is a possible future 
improvement of the present study. 
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Figure A.l: Particle spectrum of the SPSla' parameter point. 
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A. The SPSla' parameter point 



In this appendix we collect the features of the SPSla' parameter point which are relevant 



for our analysis. The spectrum of supersymmetric particles is given in Fig. |A.l| . The masses 
and decay branching fractions, generated using SPHENO fS^ and Sdecay [OS, are given in 



Tab. A.l . The consistency with low-energy constraints has been checked with SPHENO [33| 



and the cold dark matter constraint with Mi cr Omegas |35]. 
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Parameter 


Value 


Parameter 


Value 


Decay 


BR 






mh 


116.0 GeV 


xt 

Al 




53.6 % 






ruA 


424.9 GeV 


xt 




13.3 % 


mo 


70 GeV 


mr;0 


97.7 GeV 


xt 




18.5 % 


^1/2 


250 GeV 




183.9 GeV 


xt 


XiW+ 


1.3 % 


tan P 


10 


m~± 

Xi 


183.7 GeV 


if. 


A, 1 ^ 


100.0 % 


sen ri 


+ 


xt 


415.4 GeV 




' Xi^- 


92.5 % 


Ar, 


— SnO GeV 


7T) — 


1 25 3 GeV 




^ xr^^^ 


4 9 % 












^ X2^ 


2.6 /o 








172.5 GeV 


~0 
X2 


— > T-j^ T 


57.8 70 




4.69 • 10^'^ 


Tn~ 


107 9 GeV 




l/^l>^ + CC. 


15.2 % 




8.041 • 10^9 




194.9 GeV 




Z^^P^ + CC. 


11.1 % 


Ap 


1.87- 10~^ 




170.5 GeV 


x^ 




2.4 % 




0.10 




77.3 MeV 


^2 


-X?r- 


86.8% 






m) 


121.5 MeV 


^2 




4.6% 








117.4 MeV 


^2 


^ xr^r 


8.6% 



Table A.l: Details of the SPSla' parameter point. £ stands for e or /i. 



B. Comparison of EPA with the exact matrix element calculation 

Photon-induced processes at linear colliders are commonly simulated using the equivalent 
photon approximation (EPA) |^T[. The WHIZARD code allows one to use the exactly gen- 
erated matrix element (cf. Sec. 2.4) and therefore provides means to test the quality of 
EPA. In the following we will expose the main differences between EPA and the exact 
matrix element calculations (EME). As an example we take the photon- induced tau-pair 
production process with successive leptonic tau decays. 



#evl/blii 



#>n.t/l,l 




E [Gi-V] of I- E [Gi'V] of ,1+ 

Figure B.l: Lepton energy distributions from the photon- induced r pair production processes 
(left: electron, right: muon) using the equivalent photon approximation (solid, black) and the exact 
matrix element (dotted, red). 



The total cross section for photon-induced taus with their successive decay after pres- 
election at 500 GeVare as follows: 



EPA 



25.495(4) pb, a. 



EME 

7T 



21.392(70) pb, 



(B.l) 
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E [GeV] ot ,1+ 



Figure B.2: As Fig. B.l, but now plotting the difference between the EME and the rescaled 
EPA, iV,™E _ jYrcs_ rpj^g error from the exact result is also given (dotted, red). 



Process 


presel 


no p± 


no E 


no 9 


no A(j) 


all cuts 


exact 


787927 


22685 


10 


14 


867 


10 


EPA 


767322 


27974 






1546 





Table B.l: Comparison of cut efficiencies as in Table || for both EPA and exactly generated 
photon-induced r's, normalized to 1 million events before preselection cuts. 



so the EPA overestimates the number of generated taus. Figure BJ. shows the difference 
between lepton energy spectra generated with EPA and EME for electrons (left) and muons 
(right). The shapes of the distributions are roughly the same, so one may be tempted to 
account for the larger EPA distribution by a simple rescaling factor K = a^^^/a^^^, 
i.e. taking Nf^^ = K Nf^^ in the i-th bin. Figure B.2 shows the difference between 
EME and such rescaled EPA distributions demonstrating that the rescaled EPA clearly 
underestimates the number of low energy leptons. For investigations that rely on the 
lepton spectra at low energies this feature becomes crucial. In our analysis, some of the 
cuts were tailored to remove the low-energy leptons from photon-induced background, in 
particular the p_\_ and cuts. Table |B.1| shows the response of a sample of 10^ EPA 
and EME events each to the cuts. Employing the EPA instead of EME calculations would 
lead to an artificial reduction of the SM background by roughly 13%; for an integrated 
luminosity of 1 ab^^ this corresponds to a 5a" deviation from the expected SM background. 
It is therefore crucial to use the exact matrix elements in calculations of photon-induced 
processes. 
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